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Summary 

A  new  class  of  parametric  model?  for  the  spectrum  of  a 
scalar  time  series  is  proposed,  in  which  the  logarithm  of  the 
spectral  density  function  is  represented  by  a  finite  Fourier 
series.  Two  alternative  parameter  estimation  procedures  are 
described,  and  the  use  of  a  fitted  model  to  provide  forecasts  of 
future  values  is  discussed.  The  model  has  been  compared  with  the 
moi.e  conventional  autoregressive/moving-average  model,  and  the 
results  of  their  comparison  are  given. 
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1.  Introduction 


Although  there  has  been  much  discussion  of  the  fitting  of 
parametric  models  to  time  series,  essentially  only  one  such  model 
or  class  of  models,  has  been  considered.  This  is  the  auto- 
regressive/moving-average  model,  which  is  discussed  by  Box  and 
Jenkins  (1970,  Chapter  3),  for  instance.  In  this  model,  it  is  , 
assumed  that  the  observed  series  {X^}  is  related  to  an  unobserv 
ed  white  noise  series  (€t>  by 

i 

Xt  =  f^A-r  +  et  “  ?ieret-r  *  (1) 


The  spectral  density  function  of  the  series  is  then 


g(ui>  =  Yi 


r=l 


i-  I  v1™ 

r=l  r 


2 

where  t  is  the  variance  of  et» 

The  model  discussed  in  this  paper  was  motivated  by  the 
observation  that  the  logarithm  of  an  estimated  spectral  density 
function  is  often  found  to  be  a  fairly  well-behaved  function,  and 
could  thus  be  approximated  by  a  truncated  Fourier  series.  The 
form  which  has  been  chosen  for  this  model  is 


g(w)  =  irr  exp(2  l  0  cos  r w)  .  (2) 

r=l  r 

The  evenness  of  g  restricts  the  exponent  to  be  a  cosine  series, 
which  is  written  in  this  form  for  later  convenience. 


It  is,  of  course,  possible  to  approximate  any  well- 
behaved  function  by  a  truncated  Fourier  series.  However, 
the  model  (2)  would  be  of  use  only  if  p  were  small.  This 
raises  the  question  of  whether  such  spectra  arise  in  practice. 
An  answer  to  this  question  will  evidently  be  found  only  in 
extensive  application  of  the  model.  The  model  (1),  used  with 
small  values  of  p  and  q  ,  has  been  remarkably  successful 
at  fitting  practical  data.  In  many  cases,  a  similarly  good  fit 
could  be  exnected  with  the  model  (2). 

In  section  2  we  discuss  maximum  likelihood  estimation  of  the 
parameters  of  the  model,  and  the  asymptotic  distribution  of 
the  estimates ,  while  in  section  3  we  consider  some  inefficient 
estimates.  The  construction  of  predictions  of  future  values 
of  the  series  is  described  in  section  4.  The  model  was  fitted 
to  a  number  of  time  series  discussed  by  Box  and  Jenkins  (1970, 
p.  575),  and  the  results  are  contained  in  section  5.  Finally, 
in  section  6  we  examine  briefly  the  effect  of  one  type  of 
departure  from  our  assumptions. 
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2.  Maximum  Likelihood  Parameter  Estimation 

Walker  (1964)  discusses  the  estimation  of  a  general  model  of 
the  form 

2 

g(u>)  s  ^jh(u;0) »  * 

where 

h(«;0)  =  |1  +  l  a  (0)eirwj 2 
r=l  r 

Now  this  condition  is  equivalent  to 


/*  log{h(w;0)  }du  =  0, 

-tr 

for  all  values  of  0  ,  and  hence  the  model  (2)  falls  into  this 
class . 

Walker  shows  that,  if  the  series  is  Gaussian,  then  the  log- 
likelihood  of  the  parameters,  given  observed  data  x^,...,x  ,  is 
approximately 


-n  log(  2n)  -  £log  t2  -  , 

-TT 


where 


i  I  n 

I(u)  =  (27m)“  l  x.e 
1  r=l  z 


.itwl 


the  periodogram  of  x^,...,xn  .  Thus  0  may  be  estimated  by 
minimising 


fn  I(w) 

ETuTST 

f  •* 

-n 


dto  , 


and  r  may  be  estimated  by  the  minimised  value. 


<3> 
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Walker  also  shows  the  resulting  estimates,  0  and  f  ,  are 

consistent,  and  furthermore  that  the  distribution  of  /n(0  -  8„ ) 

converges  to  N{0,  WO^)"1}  where 

%  %  %u 

W(C)  5  f  log  feCw;0)>{T|  log  h(o>;8)  3  dco  • 

*  ^  - 
—  TT 


For  the  model  (2),  it  is  easily  verified  that  W(0)  is 
identically  equal  to  I  ,  which  is  one  reason  for  the  use  of 
2  cos  rw  in  that  model. 

With  other  models,  direct  minimisation  of  (3)  is  not  the  usual 
procedure  for  obtaining  the  estimates  0  and  t2  .  However,  for 
the  model  (2),  the  result  that  W(0)  =  I  makes  such  a  procedure 
simple.  For  a  Newton-Raphson  minimisation  gives  rise  to  an 
iterative  scheme  in  which 


A  .  A  1  A 

e  -  v(0  )  v( 0  ) 
-n  —  -n  n 


* 


A  A 

where  *v(8n>  and  ViS^)  are  the  vector  and  matrix,  respectively, 

o*  first  and  second  derivatives  of  (3),  evaluated  in  each  case 

at  0  .  Now  V(§  )  can  be  approximated  by  V(0n),  which  can 
~n  ~  ~n  -  -u 

further  be  approximated  by  its  expectation,  since  its  variance 
is  readily  seen  to  be  0(n“*).  But  E{V(0q)}  may  itself  be 
approximated  by  2t2W(0q)  =  2t2I  .  Replacing  t2  by  its 
n'th  iterated  estimate,  f*,  which  is  of  course  just  (3)  evaluated 
at  0n,  we  obtain  a  modified  iterative  solution,  in  which 


0 

~n 


v  ( 9  )  . 
~  n 


s 


Note  that 


ll 


I  ; 

-IT 


2ttI(m) 

*h(w;§  ) 
n  n 


c(w)dcL) 


where  c(w)'  =  (cos  w, . . . jcos^pw) .  Now  the  first  factor  in  the 
integrand  is  just  I(a>)/gn(u>) ,  where  gn(u>)  is  the  n'th 
iterated  estimate  of  the  spectral  density  function  g(o>)  .  Thus 
the  vector  of  modifications  contains  the  Fourier  coefficients  of 
the  normalised  periodogram,  and  in  particular  these  vanish  at  the 

A 

minimising  values  0  . 
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3.  Inefficient  Parameter  Estimates 

Some  consistent  but  inefficient  estimates  of  the  parameters 

of  our  model  are  suggested  by  the  following  observation.  Let 

*  2irj/n,  and  m  *  C(n-l)/2].  Then  the  distribution  of  the 

periodogram  ordinates  I(w,  ),...,I(u.)  ,  may  be  approximated  by 

i  n 

that  of  m‘  independent  exponential  random  variables  with  means 


g(w.)  =  ^  exp(2  I  0r  cos  ru.)  . 

r*  1  ^ 

Suppose  that  Z^,...,Zm  have  precisely  this  distribution,  and 
let  Cj  =  log  Zj  ,  j  =  l,...,m.  Then  the  Cj  ere  uncorrelated. 


E(Cj)  *  i|»(l)  ♦  log^j-  +2  l  ®rco8  rwj  , 


var(c^)  =  ip •  ( 1)  . 

Here  is  the  di-gamma  function,  whence  ij>(l)  is  -y  ,  and 
y  is  Euler’s  constant,  0.57722.  Also  ip  *  Cl  >  is  ir^/6  =  1.64493. 


t‘»  =  i  i  K. 
r  n  3 


cos  r(0j  j  r  -  1  *  •  •  •  jp  9 


u(1)  =  2 it  expC-^’Cl)  +  |  ^  c.) 


Then  it  is  easily  verified  that 


E(t£15) 

varCt^1*) 

cov(t^1),0^1)) 


=  0  +  0(n-1)  , 

r 

=  ^'(1)  +  0(n-1) , 


CXn”1), 
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and  indeed  that  the  distribution  of  /n(t^  -  6)  converges  to 
N(0,f' (1)1) . 

How  if  we  define  0***  in  terms  of  log  I(«^)  ,  j  *  l,...m, 
in  the  same  way  that  t*1*  is  defined  from  ,  j  *  l,...,m, 
then  clearly  under  suitable  conditions  the  same  asymptotic 
distribution  will  be  valid.  Comparison  with  the  results  for  the 
approximate  maximum  likelihood  procedure  described  in  the  previous 
section  shows  that  the  estimates  0*^  have  efficiency  ^'(l)*1  , 
that  is,  around  61%. 

Jones  (1964)  has  discussed  the  calculation  of  quantities 
very  similar  to  e^1*,  but  recommends  that  the  periodogram  be 
smoothed  somewhat  before  taking  logarithms.  Even  in  terms  of 
Z1,...,Zm  ,  this  gives  rise  to  rather  difficult  problems  in  deter¬ 
mining  the  properties  of  the  estimates.  However,  a  similar 
procedure  would  be  as  follows.  Suppose  that  m  has  a  fairly 
small  factor  k  ,  and  m  =  kn^  .  Then  the  Z's  may  be  smoothed 
in  a  sense,  by  dividing  them  into  m^  adjacent  blocks  of  k  , 
and  summing  the  blocks.  Let 


zj 


do  . 


4: 


fc=(i-l)k+l 


<£  0  ]  •  1  y  •  •  •  fin  • 


and  ?{k)  =  log  zfk) 


Clearly  the  expectations  of  Z’s  in  the  same 
block  will  be  similar.  Suppose  they  were  identical,  i.e. 


0  cos  ru>.  = 
r  t 


J^coB  ™jk 


l  =  (j-l)k+l,...,jk 
j  =  1 , . * jm^  . 


It  would  then  be  true  that 
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E(cjk>) 

var(Cjk*) 


P  2 

s  I|»(k)  +  2  l^rcoa  rWjk  +  log^y  , 


=  tjj’  (k)  , 


j  =  .  Thus  one  may  define  new  statistics 


.(k)  .  2k 


m. 


zk  p  r  ( k  )  _  _ 

“n  jli  cj  008  rwj  » 

2 it  exp(-*'(k>  ♦  ^  l 

n  j=l  3 


r  =  l,...p, 


u 


CD  . 


m. 


,(k) 


The  asymptotic  distribution  of  these  statistics  as  n  •  for 

(k) 

fixed  k  is  such  that  the  distribution  of  /n(t  -0)  converges 

/if  \ 

to  N(0,kTj»' (k)I) .  As  before,  we  define  0  in  terms  of 

A/  %  % 

logarithms  of  sums  of  periodogram  values,  in  the  same  way  that 

/if) 

t  is  defined  from  the  Z’s.  Then  under  suitable  conditions, 

(k) 

the  distribution  of  /n(0  — 0)  will  converge  to  the  normal 

(k) 

distributions  In  particular,  the  efficiency  of  9  relative 
to  the  approximate  maximum  likelihood  solution  discussed  in 
the  previous  section  is  {k^* (k)}”1,  values  of  which  are  given  in 
Table  1. 

2 

Clearly  we  may  also  define  estimates  of  t  *  to  be  defined 

(k) 

in  terms  of  ICwj)  analogously  with  u  ,  k  *  1,2, . 

The  first  of  these  was  discussed  by  Davis  and  Jones  (1968)  and 
more  recently  by  Cleveland  (1971).  The  efficiencies  of  these  will, 
for  Gaussian  data,  follow  the  same  progression  as  in  Table  1. 
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For  non-Gaussian  data,  the  situation  is  not  so  clear.  For  instance, 
the  variance  of  an  efficient  estimate  of  t  ,  as  derived  by 
Walker  (1964)  for  a  linear  process,  depends  on  the  value  of  the 
fourth  cumulant  of  the  innovation  process.  Thus  the  argument 
we  have  used  should  not  be  expected  to  be  valid  in  that  case. 


10 


■i.  Construction  of  Predictions 

For  the  model  (1),  predictions  of  future  values  may  be 
constructed  directly  from  this  equation,  as  in  Box  and  Jenkins  (1970, 
Chapter  5).  Since  the  parameters  in  the  exponential  model  (2) 
have  no  similar  direct  interpretation,  the  construction  of 
predictions  is  based  on  more  general  theory;  see  for  instance  Hannan 
(1970,  Chapter  3).  For  a  series  (Xt>  with  spectral  density 
function  g(«),  the  central  problem  is  to  find  a  function  a(*)  of 
a  complex  variable  2  ,  satisfying 

2 

(i)  g(w)  =  yjp|a(e*u)  | 2  for  some  x2, 

(ii)  a(0)  =  1 


(iii)  a  has  neither  zeroes  nor  poles  in  the  unit  disc 
|z|  <  1. 

For  if  the  Taylor  series  expansion  of  <*(•)"*  is 


a(  z)"1  =  1  -  £  B 
s  =  l 


s 


then  the  blest  one-step  predictor  of  X^  is 


!t  \?,Bsxt-s 
8-1 


Ak  0  r)  O 

and  its  mean  squared  error,  E{(Xt  --'X^)  }  ,  is  x  . 

Predictions  may  be  made  further  into  the  future  recursively. 
For  the  best  two-step  predictor  of  X^  is  found  bv  predicting 
Xt-1  by  the  above  formula  with  t  replaced  by  t-1  ,  and  then 
substituting  the  predicted  value  into  the  original  formula.  An 
n-step  prediction  is  found  by  repeating  this  operation.  The  mean 
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squared  errors  of  these  successive  predictions  are  found  from  the 
Taylor  series  expansion  of  a(*)  .  For  if 


a(z)  =  1  + 


ct  z 


s=l  8 


> 


then  the  mean  squared  error  of  the  n-step  prediction  is 

<5  n- 1  a 

T  <1  +  I  of)  ,  n  >  1  . 
s=l  8 

For  the  effect  on  these  mean  squared  errors  of  using  fitted 

parameter  values,  see  Bloomfield  (1972). 

General  numerical  methods  are  available  for  obtaining  the 

coefficients  {<*  }  and  {0  }  .  However,  for  the  exponential 

s  s 

model  (2)  we  may  use  the  structure  of  the  model  to  obtain  these 
coefficients  analytically.  For 


a(z)  =  exp( 


P 

0 

r=l 


is  clearly  the  required  function.  Suppose  first  that  p  =  1,  so 

that  o(z)  =  expCO-^z)  .  Then  clearly  cts  =  08/s!  and 

0  =  (-1)  a_  Next  suppose  that  p  >  1  .  Then 

s  s 

p  r 

o(z)  =  n  exp(0^,z  )  . 
r=l  r 

The  Taylor  series  of  each  factor  of  this  product  is  easily  written 
down,  and  the  Taylor  series  of  the  product  is  just  the  convolution 
The  coefficients  {0^}  may  be  obtained  similarly. 

The  auto- covariance 8  of  a  series  (X^)  may  be  obtained  by 
a  similar  argument.  As  before,  we  suppose  first  that  p  =  1  . 

Then  the  spectral  density  function  is 


*i*!KM*V9t* 
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h  expte^eK  e-iw)  =  f  Is<201)(ei“)s  , 

s — — 00 

where  Ig(»)  is  a  modified  Bessel  function  (Abramowitz  and 

Stegun,  1964,  p.  374).  Hence  the  auto-covariance  of  lag  s  is 
2 

x  Is ( 2©1 )  .  When  p  >  1  ,  the  spectral  density  function  is  a 

product  of  similar  terms,  and  thus  the  autocovariances  are 

convolutions  of  sequences  of  Bessel  functions.  Notice,  however, 

that  the  convolutions  in  this  case  invdlve  infinite  summations, 

whereas  those  involved  in  obtaining  the  coefficients  {ag}  and 

(B  }  involved  only  finite  sums. 

8 

If  one  goes  back  over  the  data  which  were  used  to  fit  the 
model,  one  may  use  the  fitted  predictor  to  predict  each  value 
from  the  preceding  value,  and  hence  find  the  error  of  this  pre¬ 
diction.  Had  one  used  the  true  predictor  to  do  this,  the  errors 
would  have  been  uncorrelated  with  each  other.  Now  it  is  easily 
shown  that  the  periodogram  of  these  errors  may  be  found  approxi¬ 
mately,  up  to  a  constant  multiple,  by  dividing  the  periodogram  of 
the  data  by  the  estimated  spectral  density  function.  However,  it 
was  shown  in  section  2  that  if  the  efficient  fitting  procedure 
described  there  is  used,  then  the  first  p  Fourier  coefficients 
of  this  normalised  periodogram  vanish.  This  implies,  of  course, 
that  the  sample  autocovariances  of  the  fitted  prediction  errors 
vanish,  for  lags  1  to  p,  to  this  order  of  approximation.  An 
examination  of  the  fitting  procedure  of  section  2  shows  that  no 
other  model  with  the  same  number  of  parameters  can  achieve  this. 
Thus  of  all  models  containing  a  given  number  of  parameters,  the 
exponential  model  (2)  gives  the  whitest  residuals,  in  this  sense. 


13 

By  contrast,  when  a  p-parameter  autoregressive  model  has 
been  fitted,  the  residuals  are  uncorrelated  with  the  previous  p 
values  of  the  original  series.  Thus  in  an  autoregressive  model, 
one  detects  lack  of  fit  by  finding  autocorrelation  amongst  the 
residuals.  One  must  evidently  look  elsewhere  in  the  case  of 
fitting  the  exponential  model. 


t 
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5.  Results  of  Fitting  the  Model 

In  order  to  compare  the  exponential  model  with  the  establish¬ 
ed  autoregressive-moving  average  models ,  it  was  fitted  to  six 
series  studied  by  Box  and  Jenkins  (1970,  pp.  525  ff).  For 
several  of  these  series.  Box  and  Jenkins  first  differenced  the 
data  in  order  to  render  it  more  tractable.  In  each  case,  the 
same  degree  of  differencing  was  used,  and  then  the  same  number  of 
parameters  were  fitted.  The  results  are  presented  in  Table  2. 

The  exponential  model  was  fitted  by  the  efficient  iterative 
scheme  described  in  section  2.  Initial  values  were  taken  to  be 
the  first  of  the  inefficient  estimates  of  section  3.  The  criterion 
for  convergence  was  based  on  the  sum  of  squares  of  the  modifications 
to  the  parameters  multiplied  by  the  series  length  and  divided  by 
the  numer  of  parameters  being  fitted.  The  iteration  was  terminated 
when  this  quantity  first  fell  below  0*01,  or  in  other  words, 
when  the  root  mean  square  modification  first  fell  below  one 
tenth  of  the  common  asymptotic  standard  deviation.  The  multi- 
perameter  models  were  fitted  sequentially,  each  new  parameter 
requiring  around  four  iterations. 

The  goodness  of  fit  of  the  models  has  been  assessed  by  the 
residual  variance*  Since  this  estimates  the  one-step  prediction 
error  variance  which  would  be  achieved  by  using  the  fitted  model, 
it  provides  useful  criterion  with  a  practical  interpretation.  For 
three  of  the  series ,  the  exponential  model  provides  virtually  as 
good  a  fit  as  the  respective  best  autoregressive-moving  average 
model.  For  another  two,  however,  the  fit  is  considerably  worse. 

5 

! 

\ 

4 
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For  those  series,  series  A  and  E  ,  the  addition  of  two 
more  parameters  to  the  exponential  model  did  little  to  improve 
the  fit. 

While  these  results  contain  no  evidence  that  time  series 
exist  which  fit  the  exponential  model  better  than  the 
autoregressive-moving  average  model,  at  least  they  show  that  the 
two  models  provide  comparable  fits  for  certain  series.  Since 
there  are  p+1  distinct  autoregressive-moving  average  models, 
each  containing  p  parameters,  one  must  expect  the  exponential 
model  to  be  rather  less  flexible.  Further  experience  will  be 
required  to  determine  whether  this  lack  of  flexibility  will  be 
sufficient  to  out-weigh  the  computational  andi  statistical 
simplicity  of  the  exponential  model. 
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6.  Departure  from  Assumptions 

When  fitting  one  of  a  sequence  of  parametric  models  to  a 
timer-series,  one  is  never  more  than  trying  to  approximate  the 
structure  of  the  series,  if  indeed  it  possesses  a  stationary 
structure.  Thus  it  will  rarely  be  true  that  the  spectrum  of  the 
series  belongs  to  the  class  of  spectra  being  fitted.  Since  this 
is  one  of  the  assumptions  underlying  the  result  of  Walker  (1964) 
which  has  been  used  in  section  2,  it  is  instructive  to  examine 
what  happens  if  that  assumption  is  dropped. 

Suppose  then  that  the  spectrum  of  the  series  is  gQ(a)), 
and  that  the  class  of  models  under  consideration  is  of  the  form 

7?  h(w;0)  , 

as  in  section  2.  For  parameter  values  0  ,  the  mean  squared  error 
of  the  corresponding  one-step  predictor  is 

IT 

{go<«)/h(w;0)}dw  # 

-TT 

We  shall  regard  the  closest  member  of  the  family  as  being  defined 

o 

by  parameters  0Q ,  Tq  ,  where  0Q  minimises  this  integral  and 
tq  is  the  minimum  value.  It  may  be  shown  that  the  estimates 
defined  in  section  2  are  consistent  estimates  of  these  quantities. 

In  that  sense,  that  estimation  procedure  and  the  present  definition 
of  closeness  are  well-matched.  However,  the  asymptotic  distribution 
of  the  estimates  is  no  longer  the  same.  Whilst  asymptotic  joint 
normality  with  zero  means  still  holds,  the  asymptotic  dispersion 


where 


matrix  is  now  V(8n)“1U(9n)V(0n)“1  , 
%  %  u  %  u  %  u 


,  fit  (2h.  (o,,0)h,  (o»,0) '  H,.vO),0)l 


U(0) 
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and  h.  and  H„  denote  the  vector  and  matrix  of  first  and 

- 1  -L 

second  0-derivatives  of  h  ,  respectively. 

For  the  p-parameter  exponential  model  (2),  it  is  easily 


seen  that 


h,(o,;9)  =  2h(w;0)c(u)  , 
Ho(o,;0)  =  4h(a,;0)c(o»)c(a,) '  , 


where  c(o,)'  =  (cos  cos  pm).  Thus  in  this  case. 


gn(w) 


-if 

Y<!>  »  *  J:  c(a,)c(W)*  da,  , 

-TT 

.  _-irfg0<w)  l2cCa>)c(a>)  ’  da,  , 

y(5}  “  *  J  OTTST'  *  * 
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Then  it  is  easily  verified  that 


U<0)  -  V<0) 


since 


,rir[  gn(w).-  ]  fgn(w) 

*  J  jErUr!-!!'2,t£(“>e(“),|HTirTFr^!!(2,td“  ’ 


r» 


c(o>)c(a»)  'da,  =  ttI  . 


Since  the  integrand  is  symmetric  and  nonnegative  definite, 
it  follows  that  U<0)  -  V(0)2  and  hence  also  V(0)“1U(0)V(0)’1  - 

%  *0  *0  %  %  %  %  %  % 

are  nonnegative  definite.  '  •' 

Thus  the  asymptotic  dispersion  matrix  for  0  derived  in 
section  2  is  in  reality  a  lower  bound,  only  achieved  when  the 
spectrum  of  the  series  being  analysed  belongs  to  the  class  of 
spectra  being  fitted.  It  seems  plausible  that  the  same  should 
hold  for  an  arbitrary  model,  but  we  have  been  unable  to  extend 
the  result  beyond  the  exponential  model. 

It  follows  that  when  fitting  the  exponential  model  as  an 
approximation,  misleadingly  significant  fitted  values  will  arise 
more  frequently  than  would  otherwise  be  expected,  when  the 
true  parameter  value  is,  in  fact,  zero.  Extra  caution  will  be 
needed  to  avoid  the  inclusion  of  unnecessary  parameters. 
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Table  1.  Variance  and 

Efficiency  of 

t(k) 

<*• 

n  var  t*k>  =  k^'(k) 
r 

eff.  t(k)(%) 

t|>'  (k) 

1 

1.64493 

61 

1.64493 

2 

1.28986 

78 

0.64493 

3 

1.18479 

84 

0.39493 

4 

1.13528 

88 

0.28382 

5 

1.10660 

90 

0.22132 

6 

1.08792 

92 

0.18132 

7 

1.07478 

93 

0.15354 

8 

1.06504 

94 

0.13313 

Series 

Table  2. 

Results  of 

fitting  the  model 

Variance 

Degree  of 
differencing 

Number  of 
parameters 

Residual 

Box  and  Jenkins 

Exponential 

~a'q> 

0 

2 

0.097(2) 

0.146 

A 

1 

1 

0.101 

0.164 

B 

1 

1 

52.2 

52.2 

C 

1 

1 

0.018 

0.023 

D 

0 

1 

0.090 

0.121 

D 

1 

1 

0.096 

0.096 

E 

0 

2 

228 

297 

E 

0 

3 

218 

295 

F 

0 

2 

113 

115 

(1)  The 

labelling  is  as 

in  Box  and  Jenkins  (1970,  p.  524). 

(2)  These  values  are  taken  from  Table  7.13  of  Box  and  Jenkins  (1970, 
p.  239). 


REFERENCES 


Abramowitz,  M.  and  Stegun,  I.  A.  (19  64).  Handbook  of  Mathematical 
Functions,  National  Bureau  of  Standards  Applied  Mathematics 
Series  Number  55.  Washington:  U.S.A. 

Bloomfield,  P.  (1972).  On  the  error  of  prediction  of  a  time 
series.  Technical  Report  11,  Series  2,  Department  of 
Statistics,  Princeton  University  (February,  1972). 

Box,  G.E.P.  and  Jenkins,  G.M.  (1970).  Time  Series  Analysis, 
Forecasting  and  Control.  San  Francisco:  rfolden-Day. 

Cleveland,  W.  S.  (1971).  Fitting  time  series  models  for  predict 

tion.  Technometrics  13,  713-723. 

1  ■■  ■  ■  ■  ■  1  ^ 

Davis,  H.  Y.  and  Jones,  R.  H.  (1968).  Estimation  of  the  innova- 
v  tion  variance  of  a  stationary  time  series.  J.  Am.  Statist. 
Ass.  63,  141-9. 

Hannan,  E.  J.  (1970).  Multiple  Time  Series.  New  York:  Wiley. 

Walker,  A.  M.  (1964).  Asymptotic  properties  of  least  squares 
estimates  of  parameters  of  the  spectrum  of  a  stationary*' 
nondeterministic  time  series.  J.  Austral.  Math.  Soc.  4, 
363-84.  -  ~ 


